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We compare the propagation of iron and proton nuclei above 10^^ eV in a structured Universe 
with source and magnetic field distributions obtained from a large scale structure simulation and 
source densities ~ 10~^Mpc~'^. All relevant cosmic ray interactions are taken into account, including 
photo-disintegration and propagation of secondary products. Iron injection predicts spectral shapes 
different from proton injection which disagree with existing data below ~ 30 EeV. Injection of 
light nuclei or protons must therefore contribute at these energies. However, at higher energies, 
existing data are consistent with injection of pure iron with spectral indices between ~ 2 and ~ 2.4. 
This allows a significant recovery of the spectrum above ~ 100 EeV, especially in the case of large 
defiections. Significant auto-correlation and anisotropy, and considerable cosmic variance are also 
predicted in this energy range. The mean atomic mass A fluctuates considerably between difi'erent 
scenarios. At energies below 60 EeV, if the observed A > 35, magnetic fields must have a negligible 
effect on propagation. At the highest energies the observed fiux will be dominated by only a few 
sources whose location may be determined by next generation experiments to within 10 — 20° even 
if extra-galactic magnetic fields are important. 

PACS numbers: 98.70.Sa, 13.85.Tp, 98.65.Dx, 98.54.Cm 



I. INTRODUCTION 

Ultra-high energy cosmic rays (UHECRs) are particles 
of energy > 10^^ eV, which give rise to spectacular air 
showers spreading over square miles when reaching the 
Earth's lower atmosphere, but whose properties are at 
the moment poorly known because of low fluxes scaling 
roughly as E^'^. A recent review on the experimental 
and theoretical aspects of this topic can be found, for 
example, in Ref. P]. 

As the Auger Observatory is still in a phase of de- 
ployment and early analysis, most of the current data 
come from two experiments relying on different tech- 
niques, the High Resolution (HiRes) "Fly's Eye" fluo- 
rescence telescopes and the AGASA ground array. The 
first technique consists in observing the development of 
the longitudinal extent of the air shower through the at- 
mospheric fluorescence yield, whereas arrays reconstruct 
the lateral development of the shower, by detecting sec- 
ondary products reaching the ground level. 

Because of large experimental systematics in energy 
determination for both methods, and because of low 
statistics at the highest energies, the flux of UHECRs 
is still poorly known. Whereas AGASA does not see 
any cut-off in the energy spectrum even at 10^" eV 3J, 
recent results of HiRes indicate a cut-off above around 
1019-8 eV a. 

The nature of the primary particles is even less clear. 
AGASA sets an upper limit on the photon fraction of 
< 28% at 10^9 eV |3, but it is experimentally difficult 
to discriminate between light and heavy nuclei: vari- 
ous methods give results which are largely dependent 
on hadronic models and on experimental uncertainties. 



leading to discrepancies in the results 0. 

Concerning the arrival directions, all the data accu- 
mulated until now are roughly consistent with isotropy, 
at least on large angular scales. However, particularly 
at the highest energies, the lack of statistics still allows 
a substantial large-scale anisotropy of the UHECR sky. 
On small scales (a few degrees), AGASA has detected 
a clustering of events [3 but its statistical significance 
remains questionable [g]. Furthermore, the HiRes exper- 
iment didn't confirm this feature 0, but the integrated 
aperture of HiRes stereo data is still lower than that of 
AGASA. 

The Pierre Auger Observatory which is currently under 
construction will within a few years accumulate a statis- 
tics on UHECRs at energies above 10^^ eV which will be 
orders of magnitude above previous experiments. The 
resulting spectrum, arrival direction map and perhaps 
composition should therefore put serious constraints on 
UHECR origin and propagation models. This provides 
a motivation for investigating a wide variety of UHECR 
scenarios consistent with current data. 

Various exotic models for the UHECR origin have been 
proposed ^3 j particularly in light of the possible absence 
of the Greisen-Zatsepin-Kuzmin (GZK) "cut-off" [ij as 
claimed by AGASA. However, the experimental situa- 
tion about the spectrum remains unsettled and classical 
"astrophysical" scenarios remain completely plausible. 
These models are based on particle acceleration at shocks 
in powerful extragalactic astrophysical objects, ranging 
from compact objects such as 7— ray bursts (GRBs) to 
large-scale radio lobes of active galactic nuclei (AGNs), 
see for example Ref. They have the obvious ad- 

vantage to be exclusively based on known physics and 
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astrophysical objects. They should therefore be tested 
extensively before a possible rejection in light of more 
speculative scenarios. 

Much work has already been carried out in the frame- 
work of astrophysical scenarios. Various models for 
sources and particle propagation predict different observ- 
ables such as the spectrum or anisotropics of UHECRs 
observed on Earth. A major feature of UHECR prop- 
agation models is the strength and extension of extra- 
galactic magnetic fields (EGMF) which can give rise to 
important deflections depending on particle charge and 
energy. The EGMF are mostly unknown at present 
and different models for these fields can lead to differ- 
ent predictions, as the comparison between Refs. [l3.IT5l | 
and Ref. shows. In Ref. [Til, the authors use mag- 
netic fields derived from a cosmological large scale struc- 
ture (LSS) simulation with magnetic fields generated at 
the shocks that form during LSS formation, whereas in 
Ref. ^3 and Ref. (1^] fields of "primordial" origin have 
been considered. While the different models for initial 
magnetic seed fields produce different large scale mag- 
netic field distributions and, therefore, lead to different 
predictions for UHECR deflection, there is still a sig- 
nificant discrepancy between Ref. [Ij, [I3 and Ref. [Tg , 
hinting that other technical reasons may play a role here. 

The goal of this article is to extend previous work on 
proton propagation in Ref. Q to the case of iron sources. 
This allows to fill a gap in the range of plausible mod- 
els before Auger data will discriminate between various 
possibilities. The Hillas criterion shows qualitatively 
that heavy nuclei can be accelerated to higher energies 
than protons because their gyroradii in the magnetized 
accelerator are smaller. Injection of heavy nuclei can 
therefore plausibly contribute to the observed UHECR 
fiux at the highest energies. 

Heavy nuclei propagation has already been consid- 
ered previously [ig- Discrete sources at given distances 
were studied with and without unstructured magnetic 
fields A model with discrete and continuous source 
distributions and magnetic fields with a uniform Kol- 
mogorov distribution was developed in Ref. 20]. In 
Ref. 21] , heavy nuclei propagation was studied for mag- 
netized individual sources. 

In the present paper, we consider the following sce- 
nario, largely inspired by the "best-fit model" which was 
presented in Ref. 0]: The sources are distributed ac- 
cording to the baryon density obtained from a LSS sim- 
ulation of a typical local universe Their density. 
Us — 2.4 X 10~^Mpc~'^, corresponds to average distances 
between sources of ~ 40Mpc, comparable to typical 
UHECR interaction lengths at GZK energies. There- 
fore, neither the continuous source approximation nor 
the "universal spectrum" discussed in Ref. 2^ is ap- 
plicable to this case. Source densities of this order of 
magnitude are motivated by i) comparable densities of 
candidates for powerful accelerators such as AGNs and 
a) the fact that they appear necessary to explain the 
small-scale clustering observed by AGASA | m l2i . l25j]. 



The observer is located in a void next to a supercluster 
of matter, to mimic our local extragalactic environment. 
The magnetic fields there are weak, although an experi- 
mental estimate of such a quantity in our immediate sur- 
roundings (i.e. within a few Mpc) is lacking. However, 
magnetic fields of the order of a /iG in galaxy clusters are 
present and sufficient to significantly deflect high energy 
charged particles. 

We do not take into account galactic magnetic fields in 
the work presented here, although they might play a sig- 
nificant role in the degradation of anisotropy signals for 
iron even at super-GZK energies. Previous studies, e.g. 
Ref. [2^ , show that galactic and halo fields spread arrival 
directions by an angle depending on energy, composition 
and arrival direction relative to the galactic center. 

The configuration studied, namely rare sources, struc- 
tured magnetic fields, and the presence of nuclei at 
injection, complicates the simulations. The numerical 
techniques and difficulties are presented in Sect. II. In 
Sect. HI we turn to the results obtained with proton in- 
jection, and in Sect. IV we study iron injection in the case 
of negligible deflection. Both Sects. HI and IV describe 
useful reference cases. The scenario with iron injection 
including the EGMF obtained from the LSS simulation of 
Ref. Il^] is developed and interpreted in Sect. V. Finally, 
we conclude in Sect. VI. 



II. SIMULATIONS AND METHODS 

The numerical framework for our simulations is similar 
to the one described in some details in Ref. ^^l- There- 
fore we will here only remind the reader of the major 
features of the technique, and present some specific re- 
marks concerning problems raised by the study of heavy 
nuclei propagation. 

A. Magnetic fields and propagation 

The cosmic environment, i.e. the magnetic field as 
well as the baryonic density which defines the source dis- 
tribution are the same as described in Ref. 113. They 
have been computed according to a simulation of large 
scale structure formation. An extensive discussion of the 
numerical mode ling of the magnetic field was presented 
already in Ref. [Tj, [13 . Here it suffices to say that the 
magnetic fields are generated at cosmic shocks according 
to the Biermann battery mechanism and are then renor- 
malized at the end of the simulation so that for a Coma- 
like cluster the average magnetic field in the core would 
be of order of a /iG. An alternative mechanism to gener- 
ate magnetic fields at shocks is provided by the Weibel 
instabilitv |27|. According to recent investigations of this 
process |28l |29j, magnetic fields with strength amount- 
ing to a sizable fraction of the thermal energy can be 
produced at cosmic shocks on very short time scales (of 
order of the inverse of the electron plasma frequency) . In 
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any case, an important feature of our LSS simulation is 
the existence of ~ 0.01 — 1/iG scale EGMF extending over 
scales of several Mpc in and around galaxy clusters. In 
contrast, EGMF in the voids are ^ 10~^ /iG, negligible 
for UHECR deflection. Importantly, the resulting EGMF 
is consistent with statistics of existing Faraday Rotation 
Measures with lines of sight through filaments, despite 
the fact that the magnetic field strength can be close to 
equipartition value with the total energy of the gas [s^ ■ 
We therefore have a strongly structured EGMF, different 
from the idealized case of a field with a Kolmogorov spec- 
trum with spatially constant parameters used in Ref. [l9j . 

Our EGMF also differ from those resulting from uni- 
form initial seeds, in that the latter appear to be more 
concentrated in the core of collapsed structures ■ 
However, while a more concentrated field produces less 
deflection of UHECRs, at least the results from the model 
in Ref. 13] indicate that such as difference is not dra- 
matic. Different numerical models predict different re- 
sults though. Unfortunately, the behavior of the mag- 
netic field with distance from the cluster center is not 
known. The current data indicate that /iG strong mag- 
netic fields extend out to at least ~ 1 Mpc and pos- 
sibly to larger distances [s^l ■ At distances above 1 Mpc 
from a cluster core, however, probing the magnetic fields 
becomes extremely difficult because the Faraday Rota- 
tion Measure loses sensitivity in low density regions. Fur- 
thermore, the intracluster magnetic field topology is also 
poorly known, although the situation will likely improve 
in the future. 

The size of the LSS box used in our simulations is 
~ 75 Mpc, with a grid of 5123 comoving cells. The ob- 
server, modeled as a sphere of 1.5 Mpc radius, is placed 
at the border of a void, not far from a massive and mag- 
netized structure which mimics the Virgo cluster. To 
enhance CPU efficiency and allow at the same time to 
have particles reaching the observer from regions further 
than 75 Mpc, we use periodic boundary conditions in the 
simulation to duplicate the allowed propagation region of 
UHECRs. 

To take into account the "cosmic variance" which 
arises because of various possible source locations rela- 
tive to the observer, we simulate different realizations, 
for which the position of 10 sources are chosen at ran- 
dom within the box, with probability proportional to 
the baryon density. Iron nuclei or protons are injected 
isotropically from each of these random sources in the 
simulation. 

Nuclei interactions are treated as described in Ref. 
We take into account photo-disintegration, pion produc- 
tion and pair production on the low energy photon back- 
grounds. Defiection is computed by solving the Lorentz 
equation of motion: For this purpose we use a Bulirsch- 
Stoer integrator with adaptive stepsize. 

In order to measure properly the all-particle spectrum 
at the observer position, we need to keep track of ev- 
ery secondary particle during the propagation. One nu- 
cleus emitted by a source can therefore generate up to 56 



nucleons after a propagation time depending mostly on 
primary energy. When recording an "event" , no distinc- 
tion is made between primary and secondary particles, 
but the properties of the particle, namely its charge and 
mass, are recorded. 

After propagation, the analysis of simulated events is 
performed for each simulated data set, roughly in the 
same way as in Ref. |l4j : Taking into account fiuctu- 
ations in the spectra and intensities of the sources, we 
build observed energy spectra, average composition plots, 
deflection histograms, typical sky maps, auto-correlation 
functions and full-sky angular power spectra (partial sky 
coverage being in theory invertible as shown in Ref. (33|). 
In the present study we neglect any finite experimental 
resolution in energy and arrival directions. The "cosmic 
variance" associated to fluctuations in source positions 
and properties within the scenario is also computed for 
these observables. It is in general deflned as the one- 
sided median deviation from the average of the consid- 
ered quantity. The fluctuations of the source properties 
are chosen as in Ref. Each source is characterized 

by a spectral index and a luminosity Qi with 

^ oc Qr2-2 for 1 <Q. < 100, (1) 

— ^ — const, for (a) — 0.1 < cti < {a) + 0.1 . 
dai 

Here, {a) will be chosen to best fit the observed spec- 
trum. Finally, we assume that all sources accelerate to a 
common maximal energy i^max for which we will choose 
different values. 



B. Numerical Difficulties 

Our simulations allow detailed simultaneous predic- 
tions of various observables in specific scenarios. The 
drawback is of course a large CPU time consumption: 
The ratio of simulated trajectories over recorded events 
is large, of order 1000 or more depending on the consid- 
ered scenario. In the case of nuclei propagation, and for 
the EGMF we use, the recorded event yield is particu- 
larly low because of large deflections and the necessity to 
follow secondary nuclei. 

Unfortunately, backtracing the particles from the ob- 
server does not allow to predict observables for the model, 
since we do not know in advance the spectrum, composi- 
tion, and effective source distribution of observed events, 
not to mention the stochastic nature of cosmic ray inter- 
actions. 

The consequence of CPU limitations is the limited 
statistics of simulated events. As we want to simulate 
many source realizations for different scenarios (typi- 
cally around 100 source location realizations for each sce- 
nario), we restricted ourselves to 10^ events per realiza- 
tion, above either 10 EeV or 40 EeV. This is sufficient as 
a flrst step, since our goal is to explore the widest range 
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of possible scenarios rather than to focus on one specific 
model, which we might choose later in the light of Auger 
results. In the future, we plan to carry out more de- 
tailed simulations using a parallel version of the cosmic 
ray propagation code currently under development. 



1. Event Reweighting 

For CPU efficiency, we inject particles at the sources 
with a uniform distribution in the logarithm of energy. 
To predict obsorvablos with fluctuating source luminosi- 
ties Qi and spectral indexes we reweight each simu- 
lated trajectory with a factor that depends on the source 
power and the injection spectrum. This reduces the ef- 
fective number of events for statistical quantities such as 
histograms and the anisotropy observables. Reweighting 
can also bias the event maps by causing spots to appear 
on the maps. 

We show here as an example how applying weights 
to simulated trajectories reduces the sensitivity of the 
angular power spectrum to large-scale anisotropics. 

• Absence of weights. Given N arrival directions Hi 
distributed on the sphere, we define our estima- 
tor of the Ci by aim = N~'^J2iLi'^em{ni) and 
Ce = (2^+ l)-iEL=-£|a^m|': where F,„(n,) 
arc the usual spherical harmonics. For the null hy- 
pothesis of full isotropy, the mean value for this 
estimator becomes 



^ ' m i,j 



As the arrival directions are independent, only the 
terms i = j contribute, for which {\Y£m{ni)\^) = 
(47r)~^/ dn\Yimin)\'^ = (47r)"\ and therefore, one 
finds the well-known result 
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(2) 



• Effect of weights: Let \is now assign weights uji to 
the events, with a distribution p{u)) with mean n 
and standard deviation a. The first average to con- 
sider is over arrival directions. The computation is 
the same as before, but since the definition of the 

estimator is now a^m = (j2j^j^ J2iL0iYirn{ni), 

it follows that (Q) = {4tt)-^ {J^ojf) ■ 
The next step is to average over the weight distri- 
bution, which is independent of arrival direction 
distribution. We place ourselves in the limit of 
large N, therefore replacing the sums over by 
integrals. This leads immediately to 



where N,, 



N 



1 



2 • 



(3) 



This last formula shows that if the weight distribution 
is broad, i.e. a/fj, ^ 1, the effective number of arrival 
directions becomes small as the only events which will 
be "counted" arc the ones with the largest weights. This 
results in an increase of the bias in the power spectrum 
estimate in the isotropic case, and the sensitivity to pos- 
sible small anisotropics in the model is reduced. 



2. Effect of Finite-Sized Observer and Simulation Box 

The ideal observer in the method developed here 
should be point-like, but, once again for CPU reasons, 
we modeled the observer as a sphere of radius 1.5 Mpc. 
There is in general no problem with this procedure, as 
we expect the UHECR properties such as density, spec- 
trum, composition, and anisotropies to be roughly the 
same within this volume. This is because the sources are 
extragalactic and typically located at distances of tens 
to hundreds of Mpc, and the magnetic field around the 
observer is very weak, ^^10 pG. However, in some real- 
izations for the source locations, one or a few sources can 
be located within a few Mpc from the observer, and the 
finite size of the observer "spreads" the image of such 
sources over a few degrees. This smoothes out the auto- 
correlation function in the first few bins of small angular 
off-set. For the case of propagation in a regime with 
considerable deflection, this effect is, however, of little 
importance. 

The finite size of the simulation box can also lead to 
spurious effects observable on simulated sky maps. For 
events at the lowest energies, E 10 EeV, when particles 
travel over large distances and in the absence of deflec- 
tions by magnetic fields, an excess of events is observed 
in the direction of the simulation box corners because of 
the rectangular geometry of the box, and the periodic 
repetition of sources can be observed in the sky maps of 
arrival directions within a given source realization. How- 
ever, for the cases of interest in our simulations, energies 
E > 40 EeV and/or non- negligible magnetic fields, these 
effects disappear completely. 



3. Specificities of Heavy Nuclei Propagation 

There are a number of difficulties which arise when 
studying heavy nuclei and their secondaries in the frame- 
work of our method, in particular concerning the analysis 
of angular distributions for simulated events. 

• Propagation of nuclei takes more CPU time 
than protons because the deflections are larger, 

which leads the Bulirsch-Stoer integrator to choose 
smaller stepsizes for integrating the trajectories. 
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FIG. 1: Histogram of injection energies of events detected by 
the observer in a scenario with iron injection. The injection 
spectrum is flat in the logarithm from 10 EeV to 10 ZeV, but 
evidently this is not the case for the detected events. This 
effect is explained in the text below. 
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FIG. 2: Auto-correlation of events with energies above 80 
EeV for a realization of iron injection with magnetic flelds, 
without (black) and with (red) correction for nuclear showers 
arising from the disintegration of a single iron near the ob- 
server. M{9) — 1 (dashed line) corresponds to an isotropic 
distribution. 



• For heavy nuclei injection, a consequence of photo- 
disintegration processes is that particles with a 
given observed energy have an extremely large 
range of injection energies. Fig. ^ shows the dis- 
tribution of injection energies for recorded events. 
Two populations appear clearly: 1) At low injec- 
tion energies, the flat distribution of iron originat- 
ing from nearby sources which is the same as in the 
case of proton injection; 2) Nuclei emitted at the 
highest energies which are detected as light nuclei 
after having undergone photo-disintegration. Due 



to the approximate conservation of the Lorentz fac- 
tor in interactions, this population appears at in- 
jection energies E > A x Emim which corresponds 
to \og{E/EeV) > 2.75, as seen in Fig. [T] Note 
that the dip between these two regions is an artifact 
due to the fact that iron at those energies undergo 
frequent photo-disintegrations but their secondary 
protons have energies below -Emin and, therefore, 
are not detected. This affects the reconstruction of 
the injection spectrum in that region. 

As events are reweighted according to their injec- 
tion energies, the situation depicted above leads to 
a broad weight distribution. Therefore, the prob- 
lems raised previously are particularly acute for 
iron injection. 

• Although they are necessary to evaluate the spec- 
trum and composition, "showers" of nuclei which 
are generated from a single iron near the observer 
create fake signals in the auto-correlation. A crude 
but efficient approach to this problem is to use only 
the first particle of each shower (composed of 56 
particles at most) in the anisotropy analysis, or 
equivalently to assign to each particle a weight in- 
versely proportional to the size of the shower it be- 
longs to. An example is shown in Fig. |21 where 
the auto-correlation is computed with and without 
reweighting of showers, using the estimator 



LJiUjj if Oij in the bin of 9 







otherwise 



(4) 



Here, S{9) is the solid angle size of the correspond- 
ing bin, and ui are weights as before. The nor- 
malization factor C = r^e/ WiWj, with llg de- 
noting the solid angle of the sky region where the 
experiment has non-vanishing exposure, is chosen 
such that an isotropic distribution corresponds to 

N[e) = 1. 

We point out that in terms of CPU time consumption, 
the scenarios explored in this paper are perhaps the most 
challenging for obtaining sufficient statistics. This is the 
result of combining heavy nuclei propagation with con- 
siderable magnetic fields and a discrete distribution of 
rare sources. Indeed, these kinds of scenarios have not 
been studied yet in great detail because they are compu- 
tationally time-consuming. However, it is important to 
perform simulations and to predict observables for those 
scenarios as they allow to scan a large range of UHECR 
models. 



III. WARMING UP: PROTON INJECTION 

To have a reference and compare with iron simulations, 
we extensively simulated proton emission and propaga- 
tion with parameters similar to the "most-favored" sce- 
nario number 6 of Ref. In the framework of this 
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TABLE I: List of simulations carried out. Each simulation 
corresponds to 10* recorded "events" . To increase the statis- 
tics at high energies, some scenarios were simulated twice: 
once with Emin = 10 EeV, and once with iSmin = 40 EeV. 

scenario, for a source density of 2.4 x lO^^Mpc^"^, we 
chose a distribution of sources following the baryon den- 
sity, with and without EGMF. The maximum injection 
energy was set to 1 ZeV= 10^^ eV, and all events above 
^-min = 10 EeV are taken into account. The injection 
spectrum in this energy range was chosen oc where 
the average value of a is fitted to the data and typically 
varies between ~ 2 and ~ 2.4. 

The main conclusions from these simulations are: 

• The predicted energy spectra are shown in 
Figs. 121 and ^ As expected, in this astrophysical 
scenario the GZK feature appears at energies above 
4 X 10^^ eV. Since the sources are distributed 
at various distances from the observer, the feature 
does not appear as a sharp cut-off, but rather as a 
steepening of the spectral index. Let us remark 
that we are not able to study the ankle within 
these simulations which are limited so far to ener- 
gies above 10 EeV. A proton average injection spec- 
tral index (a) = 2.4 fits better than (a) — 2.0 to 
current AGASA/HiRes data at 10 < < 40 EeV, 
see Fig. 0] 

A major aspect of the obtained spectra is the pres- 
ence of a large cosmic variance, which has differ- 
ent origins below and above the GZK feature: Be- 
low the GZK feature, this variance is mostly due 
to source luminosity fluctuations, since it almost 
disappears if sources have identical properties, see 
Fig. 21 Above the GZK feature, fluctuations of the 
source locations also significantly contribute to cos- 
mic variance: The spectrum at post-GZK energies 
depends substantially on the presence of a few pow- 
erful sources close to the observer. This is because 
our source density is low: No continuous approxi- 
mation for the source distribution could reproduce 
such results. 

Finally, the comparison of the upper and lower pan- 
els in Fig. demonstrates the effect of EGMF on 
the spectrum: The attenuation of spectra at high 
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FIG. 3: Energy spectra predicted in the case of proton in- 
jection at the sources, without (top panel) and with (bot- 
tom panel) EGMF. The average (solid line) and cosmic vari- 
ance (shaded band) result from various realizations of the 
source properties and locations. The source positions are as- 
sumed to have a statistical distribution proportional to the 
baryon density. For each source location, 50 realizations for 
source intensity and spectral index are drawn according to 
the distributions of Eq. Q . The mean source spectral index 
is (a) = 2.4. To guide the eye we also show the spectra mea- 
sured by AGASA W\ and HiRes Q. The normalization of our 
curves is obtained by adjusting the average spectrum to low 
energy data. 

energies is slightly more pronounced in the case of 
EGMF, the mean spectral slope in the energy range 
between ~ 50 EeV and ~ 200 EeV being - 6.4 in- 
stead of ~ 5.3. This is due to the increase of the 
average traveled distance in the presence of EGMF. 

• Deflections and time delays are defined here as 
the deflection and delay accumulated by UHECRs 
compared to photons. In general they are very dif- 
ficult to observe, but i) moderate deflections from 
an unambiguously identified source could be mea- 
sured, and a) time delays could be measured in 
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FIG. 4: Same as Fig. for proton injection with EGMF 
fields: In the top panel all sources have the same intensity 
Q = 1 and spectral index a = 2.4, i.e. cosmic variance is only 
due to fluctuations in the source locations. The bottom panel 
includes source property fluctuations as in Fig.|3 but with a 
mean injection spectral index (a) — 2.0 instead of 2.4. 



case a GRB or a SN explosion can be identified as 
source. Predicted deflection angle and time delay 
histograms are shown in Fig. |31 It can be seen that 
above 10 EeV the deflection distribution obtained 
is not trivial, i.e. it is not the one we would ob- 
tain if particles were in the fully diffusive regime. 
One can also remark that even at 100 EeV, typical 
deflections are still of the order of 10-40 degrees. 
This is due to the fact that EGMF are substantial, 
and sources are generally located in magnetized ar- 
eas: particles are mostly deflected in the local en- 
vironment around their sources. The time delay 
histogram shows that typical time delays are large, 
around 1 Gy at 10 EeV, due to the EGMF. Thus 
in this scenario, even at the highest energies time 
delays will be too large to be directly measurable. 

• Auto-correlation of events. Even for substan- 
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FIG. 5: Distributions of deflection angles (top panel) and time 
delays versus straight-line propagation time (bottom panel) of 
simulated trajectories for proton injection with spectral index 
a = 2.4. The distributions are cumulated over all source loca- 
tion realizations and the error bars are drawn from Poissonian 
statistics and reflect the finite number of simulated trajecto- 



tial EGMF, the fact that at high energies all events 
come from a few sources leads to strong auto- 
correlation signals. A major feature is, however, 
the high cosmic variance of this auto-correlation, 
which depends strongly on the positions and lu- 
minosities of the sources closest to the observer. 
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This can be seen in Fig. which shows two dif- 
ferent source reahzations: The first one presents a 
smooth auto-correlation extending to 30 degrees, 
whereas in the second realization, with a nearby 
source, the auto-correlation is highly peaked, with 
a gaussian-like shape of width a ^ 5 degrees. 
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FIG. 6: Examples of auto-correlation for E > 100 EeV for 
proton injection with spectral index a — 2.4 and EGMF, for 
two different source configurations. Within the same model, 
fluctuations of the distances to the nearest sources from the 
observer generate extremely different auto-correlation (and 
therefore clustering) shapes. N{8) = 1 (dashed line) corre- 
sponds to an isotropic distribution. 
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FIG. 7: Energy spectra predicted for iron injection with- 
out deflection. The average (solid line) and cosmic variance 
(shaded band) are obtained as in Fig. |3 The average source 
spectral index is (a) = 2.0, whereas the dashed lines show 
the average spectrum for (a) — 2.4. The maximum injection 
energy is set to 4 ZeV (top panel) and 10 ZeV (bottom panel). 



IV. IRON INJECTION WITHOUT MAGNETIC 
DEFLECTION 

This section illustrates general features of heavy nuclei 
propagation. We do not take into account particle de- 
flections, and, therefore, these results can be compared 
to the case of small deflections obtained in the EGMF 
scenarios of Refs. ,16, 20] . Note, however, that deflection 
can still be of the order of ^ 20° in the EGMF scenario 
ofRef. Wd\ when their results are extrapolated to nuclei. 



A. Spectrum 

In Fig. we represent the spectra obtained with an av- 
erage iron injection spectral index (a) = 2, and two dif- 
ferent maximum injection energies. The injection spectra 



which allow to fit the observed spectra tend to be some- 
what harder than in the case of protons. Clearly, how- 
ever, the spectral shape predicted at energies E ^ 30 EeV 
is not compatible with observations. At those energies, 
a lighter component, e.g. from proton injection, is nec- 
essary. Around 20-30 EeV a small "bump" is predicted 
for pure iron injection which is mostly due to secondary 
sub-GZK protons. 

At energies E ^ 30 EeV we observe a flux suppression 
due to nuclear photo-disintegration. The shape of the 
cut-off is not the same as in Fig. |31 for proton injection: 
in particular it starts at lower energies and it is flatter. 
Furthermore, there is a significant steepening around 100 
EeV, with a flattening between 10^° and 10^^ eV. The lo- 
cation of this flattening or "spectrum recovery" strongly 
depends on the source conflguration, as can be seen from 
the large cosmic variances in Fig.Q and also on the max- 
imal injection energy: Fig.[7|shows that the recovery flux 
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is statistically higher when increasing E-^ 
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B. Composition 
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FIG. 8: Predicted average atomic mass (A) as a function of 
observed particle energy corresponding to the spectra shown 
in Fig.|7|for iron injection without deflection. Mean and cos- 
mic variance are computed as in Fig. |7| with (a) — 2.0 and 
maximum injection energy of 4 ZeV (top panel) and 10 ZeV 
(bottom panel). The dashed lines show the average composi- 
tion for (q) = 2.4. 

The mean composition as a function of energy is shown 
in Fig. where the cosmic variance is due to the fluc- 
tuations in source locations and properties. At energies 
i? ^ 30 EeV, the mean mass (A) is typically higher than 
35, reflecting in fact a bimodal distribution: At these 
energies we have on the one hand protons originating 
from high-energy heavy nuclei photo-disintegration, and 
on the other hand heavy particles injected and surviving 
propagation at low energies. If we model the composition 
as a simple (II,Fe) mixture, which is done in most exper- 
imental studies at these energies, the observed (A) ~ 35 
in Fig. |H1 corresponds to ~ 40% H and ~ 60% Fe even for 



FIG. 9: Distributions of masses in the case of iron injection 
with spectral index 10 ZeV and no EGMF. 

The distributions are cumulated over various source position 
realizations and the error bars are drawn from Poissonian 
statistics reflecting the finite number of simulated trajecto- 
ries. 



pure iron injection. The dashed lines in Fig. |S1 represent 
(A) in the case of a different injection spectral index, 
(a) — 2.4 instead of 2: It is evident that the steeper 
the injection spectrum, the heavier the composition. In- 
deed, decreasing a enhances the number of high energy 
injected iron and, therefore, of low-energy proton sec- 
ondaries which implies a lighter composition at energies 
below ~ 100 EeV. 

Between ~ 40 EeV and ~ 200 EeV, the situation is less 
clear as there is a competition between two phenomena: 
On the one hand, at energies above ^ 100 EeV iron is 
efficiently photo-disintegrated, whereas below -Emax/56 
secondary protons from dissociated iron nuclei appear. 
Both effects suppress the average atomic mass (A) . This 
can lead to a bump around ~ 100 EeV provided that the 
two effects operate in different energy ranges, that is if 
-E'max/56 is well below ^ 100 EeV, as in the top panel of 
Fig.lHl In general, these effects are reduced for steeper in- 
jection spectral index a, as can be seen in Fig.|Sl because 
there are fewer secondary products to affect a larger pop- 
ulation of particles in the lower energy part of the spec- 
trum. Therefore, at these transition energies, the mean 
composition appears to be subject to large fluctuations 
depending on various parameters. 

Above ~ 200 EeV, the mean value of atomic mass (A) 
depends substantially, like the spectrum, on the source 
locations. However, in general the average (A) slightly 
increases with energy, resulting from the kinematic condi- 
tion i5max(^) = £'max(56) X A/56, as photo-disintegration 
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roughly conserves the Lorentz factor. 

We also found that at all energies considered, fluctu- 
ations of source luminosities and injection spectra make 
an insignificant contribution to the cosmic variance in the 
distribution of {A). 
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FIG. 10: Predictions for the angular power spectrum in the 
case of iron injection without any magnetic field, above 40 
(top) and 80 EeV (bottom). The dashed line, computed 
from Eq. Q, corresponds to an isotropic distribution of 10* 
weighted simulated trajectories above 40 EeV. The error bars 
represent the cosmic variance obtained by simulating various 
source positions and properties. 

In Fig. 121 we present the distributions of mass A for 
pure iron injection, in the case of a = 2 and Smax = 
1 ZeV. There are two populations of detected events, 
namely a contribution of protons on the one hand, and 
a broad distribution of heavy nuclei dominated by the 
iron group elements on the other hand. This bimodal 
distribution is also reflected in Fig. ^ and justifies to 
consider the distributions as a (H,Fe) mixture at first 
approximation, as was done above. Fig. shows a nu- 
clear cascade in which the abundance of a given nucleus 
relative to its parent nucleus is governed by the ratio of 
the photo-disintegration rate of the parent and its total 



FIG. 11: Predictions for the auto-correlation function of ob- 
served events in the case of iron injection without EGMF, 
above 40 and 80 EeV. We do not show the auto-correlation in 
the first bin and showers have been corrected for, see Fig. |5] 
The error bars represent the cosmic variance obtained by sim- 
ulating various source positions and properties. N{6) — 1 
corresponds to an isotropic distribution. 



disappearance rate. For unstable daughter nuclei, this 
ratio is in general smaller than one, which explains the 
roughly exponential behavior of this cascade. 

We note that the predicted composition discussed in 
this section is not ruled out by the still very scarce ex- 
perimental data 0. 



C. Anisotropics 

In Figs. EH and ^2 we represent the autocorrelation 
function and angular power spectrum predicted in the 
case of iron injection without any magnetic deflection. 
As for the case of spectrum and composition, fluctuations 
due to various source positions and properties are com- 
puted: They are rather large, but at energies E > AO EeV 
they are smaller than the deviation from the values pre- 
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dieted from isotropy. 

In the null hypothesis of isotropy, the auto-correlation 
M{9) is 1 by definition, and the angular power spec- 
trum is flat, at a value depending on the number 
of recorded events and on their weights, computed ac- 
cording to Eq. (PJ- The isotropic level of angular power 
spectrum is represented by the dashed lines in the plots. 

The angular power spectrum is approximatively flat as 
is expected for a distribution of point sources: The sky 
map being like a sum of Dirac 5 functions, its Fourier 
transform is flat. At 40 EeV, sources contribute up to cos- 
mological distances and therefore the Ci are compatible 
with the expected value from isotropy. At 80 EeV how- 
ever, because the number of observed sources becomes 
small, and the predicted Ci is at a significantly higher 
level than expected from isotropy. 

As expected in the absence of any magnetic field, the 
auto-correlation of events is sharply peaked at small an- 
gles, reflecting also the small source density. The auto- 
correlation signal is larger at 80 EeV than at 40 EeV, as 
the number of effectively contributing sources becomes 
smaller. However, the peak in the auto-correlation at 
small angles has a finite width because i) sources are 
concentrated in regions of large baryonic density, which 
makes them often appear next to each other, and ii) there 
is an artificial spreading of the auto-correlation over a few 
degrees due to the finite size of the observer, see Sect. II- 
B. The power spectra and auto-correlations are compa- 
rable to the analogous case with proton injection ,lj|, as 
might be expected for negligible deflection. 



V. IRON INJECTION WITH MAGNETIC 
DEFLECTION 

A. Spectrum and Composition 

The spectra obtained in presence of magnetic deflec- 
tions are presented in Fig. 1121 and can be compared to 
Fig. [71 The bump which was predicted between 20 and 
30 EeV without EGMF now appears as a smooth feature 
in a broader energy range, 10 ;$ i? < 40 EeV. There is no 
steepening of the flux at 100 EeV, comparable to proton 
injection, and contrary to the case with no EGMF. 

At higher energies, we observe a flattening of the spec- 
trum as in the absence of EGMF. The flux is higher than 
without EGMF, at least statistically: For E'max = 4 ZeV, 
we predict {JE"^) - 0.1 eVcm^^ g-i sr^i at - 500 EeV 
with EGMF, instead of ^ 0.01 eVcm"^ s^^ sr^^ in the 
absence of EGMF. 

The average atomic mass {A) as a function of energy 
is presented in Fig. E| At energies i? < 60 EeV, it can 
be seen that, depending on the average injection spec- 
tral index (a), (A) - 15 - 30 instead of {A) - 30 - 50 
in the absence of EGMF. The interpretation is that, as 
explained in Ref. [T9l |. magnetic fields increase the mean 
path length between sources and observer and the result- 
ing increase in interactions also drives up the relative pro- 
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FIG. 12: Energy spectra predicted for iron injection in the 
presence of the EGMF obtained from the LSS simulation. 
The average (solid line) and cosmic variance (shaded band) 
are obtained as in Figs. El and |7] The average source spectral 
index is (a) = 2.0. The maximum injection energy is set to 
4 ZeV (top panel) and 10 ZeV (bottom panel). The dashed 
lines represent average spectra in the case (a) = 2.4. 



portion of lighter secondaries to iron. As a consequence, 
it is very interesting to remark that a measured average 
mass of UHECRs at energies 10 < < 30 EeV larger 
than ~ 35 would imply that deflections due to EGMF are 
relatively small, independently of the nature of particles 
accelerated at the source. This can be another test for 
extragalactic magnetic fleld effects, independently from 
anisotropy studies 0|. The trend of (A) with spectral 
index (a) is the same as in the absence of EGMF. Fur- 
thermore, the shape of the mass distribution is similar 
to Fig. 1^1 except that abundances of light elements with 
A 6 tend to be increased by a factor 2-3 which corre- 
sponds to the decrease in {A) seen in comparing Fig. jSl 
with Fig. US 

At higher energies the energy dependence of (A) 
around ~ 100 EeV depends crucially on the source pa- 
rameters £'max and a, as in the absence of EGMF. The 
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FIG. 13: Predicted average atomic mass {A) as a function of 
observed particle energy corresponding to the spectra shown 
in Fig. 1121 for iron injection in the presence of the EGMF. 
Mean and cosmic variance are computed as in Fig. 1121 with 
(a) = 2.0 and maximum injection energy of 4 ZeV (top panel) 
and 10 ZeV (bottom panel). The influence of the EGMF is 
apparent from comparing with Fig.|S|for the same case with- 
out EGMF. The dashed line corresponds to a mean injection 
spectrum (a) = 2.4 instead of 2. 



EGMF can lead to additional modification of the com- 
position, as seen by comparing Figs. [TSl with ISl 



B. Anisotropies: Deflections and Sky Maps 

As expected, the histogram of deflection angles for iron 
injection, presented in Fig. shows typical deflections 
larger than for proton injection. This is particularly the 
case at the highest energies, since at those energies the 
observed particles are mostly heavy nuclei, whose local 
deflection is therefore up to 56 times larger than in the 
case of protons. As a result, typical deflections of 10° or 
more are expected even at 500 EeV. 

Sky distributions of deflection angles in the framework 
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FIG. 14: Histogram of deflection angles of observed UHECRs 
for different cuts in energy, in the case of iron injection with 
a — 2. As in Fig. |K1 the distributions are cumulated over all 
source location realizations and the error bars are Poissonian, 
reflecting the finite number of simulated trajectories. 



of this scenario are presented in Fig. 1 151 for energies above 
40 and 80 EeV. The considerable inhomogeneity of these 
maps reflects the presence of extended magnetized struc- 
tures distributed along the LSS. There is a striking dif- 
ference between these maps and the one presented in 
Ref. : Even at 80 EeV, there is still a large part of the 
sky where deflections are typically > 50° . This is due to 
at least two reasons: Apart from the larger charge of the 
UHECRs in this scenario, the EGMF in our simulation 
is more extended, as discussed previously. 

These deflections may prevent us from performing 
straightforward "UHECR astronomy" with forecoming 
experiments such as the Pierre Auger Observatories. 
However, they do not erase all the structures in the sky 
either, and in particular the expected low source density 
may allow to identify extended sources in the sky even 
within this unfavorable scenario. An example is given in 
Fig. 1161 where arrival direction maps are represented for 
a given source realization with iron injection and EGMF, 
and a statistics of ~ 150 recorded trajectories above 10^° 
eV. At 40 EeV, the sky is dominated by an isotropic cos- 
mological background, but at higher energies the back- 
ground disappears and a small number of nearby sources 
emerge. These sources are clearly visible although they 
are not point-hke. It has to be noticed that, if such a 
configuration were realized in nature, then the very few 
sources one expects to contribute at the highest energies 
can be found in the northern hemisphere, where Virgo 
is located, emphasizing the need for a northern UHECR 
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FIG. 15: Sky distribution of deflections obtained from events 
recorded by the observer in the case of iron injection with 
EGMF, accumulated over all simulated realizations. Top 
panel: Events above 40 EeV. Bottom panel: Events above 
80 EeV. The maps have been smoothed with beams of 5 de- 
grees. 

observatory. 

C. Anisotropies: Auto-Correlation and Angular 
Power Spectrum 

1. Auto-Correlation Function 

The large deflections do not prevent us from observing 
a strong auto-correlation signal at high energies, due to 
the low source density. We present in Fig. El the auto- 
correlation predicted for this scenario, computed with a 
1 degree binning. Due to the computational reasons dis- 
cussed in Sect. II, we count nuclei showers as only one 
particle and we disregard the first angular bin. It can be 
seen that above 40 EeV, the predicted auto-correlation 
signal is almost flat, whereas above 80 EeV it reaches a 
factor ~ 4 above the isotropic level, and much larger val- 
ues above 120 EeV. Fig. El also shows that the predicted 
auto-correlation signal is highly dependent on the source 
configuration, and can extend over 10 degrees or more 
due to diffusion in the EGMF. 

When comparing Fig. El to Fig. ^2 it can be seen 
that the effect of magnetic fields is to strongly reduce the 
intensity of the auto-correlation peak at small angles: at 



E > 120 EeV 




FIG. 16: Example of arrival directions of events above 40, 80, 
and 120 EeV in a particular source realization, with roughly 
150 events recorded above 100 EeV. In this realization, at 
high energy an individual source emerges progressively from 
the background. 



80 EeV, the reduction factor is ~ 10. Due to the larger 
deflection of nuclei this suppression is stronger than for 
protons. 

We remark that the galactic magnetic fields, which 
were not take into account here, will also increase this 
smoothing of auto-correlation over larger angular scales. 



2. Angular Power Speetrum 

In Figs. El and El we present examples of angular 
power spectra, and their averages and fiuctuations for 
given energy thresholds. The predictions are similar 
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FIG. 17; Auto-correlation function for events above 40, 80 
and 120 EeV in the case of iron injection with EGMF. We 
do not show the auto-correlation in the first bin and showers 
have been corrected for. The error bars represent the cosmic 
variance, computed in the same way as for spectra and compo- 
sition: Various source location configurations are considered, 
and for each one 50 source intensity and spectral indexes are 
drawn. J\f{d) = 1 corresponds to isotropy. 
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FIG. 18: Example of predicted angular power spectra for two 
source configurations for the scenario with iron injection and 
EGMF. Top panel: Events with E > 40 EeV; the correspond- 
ing event map is roughly isotropic and contains several thou- 
sand events. With reweighting, we predict {Ce) — 1.54x 10~*, 
i.e. NcH — 500 effective events, see Eq. 0; this level is 
indicated by the dashed line. Bottom panel: Events with 
E > 120 EeV. Two nearby sources appear in the corre- 
sponding map. The average (Ct) predicted from isotropy is 
7.1 X 10^'^ (dashed line), and the average {Ce) predicted by the 
simulation implies A'',, — 2.2 "effective point sources". Fur- 
thermore, Ce decreases with increasing I, which reflects the 
spatial extension of the sources. 



to those obtained in the absence of EGMF in Fig. [TUl 
The theoretical prediction for the power spectrum in the 
isotropic case is given by the expressions in Eq. (Q. As 
the large fluctuations show, the power spectrum depends 
on the energy threshold, as well as on the source config- 
uration. 

Above 40 EeV the mean angular power spectrum is al- 
most flat and consistent with predictions from isotropy; 
however some realizations exhibit some large-scale pat- 
terns visible in the low-order C^. When increasing the 
energy, the angular power spectrum remains roughly flat, 
but at levels much higher than expected from isotropy. 
Indeed, at high energies there is only a small number Ny 
of bright sources visible in any given realization, leading 
to a mean value of Ci = (AirNy)'^. The power spec- 
trum is then approximatively flat, with deviations from 
flatness caused by the observed source extension due to 
deflection in the EGMF. Above 80 EeV, the average {Ci) 
is 5 X 10"'^, leading to an average value Ny ~ 16 visible 
sources, whereas above 120 EeV we have (Ce) ~ 0.03, 
and therefore, only 2.7 sources are visible on average. 

The angular power spectrum is a useful tool to esti- 
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FIG. 19: Angular power spectra for events above 40, 80 and 
120 EeV in the case of iron injection with EGMF. The error 
bars represent the cosmic variance computed in the same way 
as in Fig. 1171 The dashed line corresponds to an isotropic 
distribution of the same number of 10^ simulated, weighted 
events above 40 EeV. 



mate the effective number of observed sources above a 
given energy, but otherwise it appears to be of quite lim- 



ited use within the scenario under consideration: At en- 
ergies below 40 EeV, deflections are too large to expect 
a large-scale signal, and at high energies the large val- 
ues predicted for Ci are only due to the finite number of 
visible sources. 



3. Detection with Finite Statistics 

When comparing predictions of auto-correlations and 
power spectra such as shown in Figs. [T7tll9l with future 
data, one has to take into account also the statistical 
error due to the finite number TVobs of the actually ob- 
served events. This has not been included in Figs. llTI - IT^ 
which thus correspond to A'obs oo. The error bars due 
to finite statistics can be assigned directly to the obser- 
vational angular power spectra or auto-correlations, as 
these errors depend on the detector aperture. 

In the same way, we have not taken into account the 
detector exposure in our anisotropy predictions, as it can 
be taken into account directly in the experimental data. 

Alternatively, the finite number of experimental events 
and non-trivial exposure functions as well as angular and 
energy resolutions can be assigned to the predictions by 
drawing mock data sets from maps constructed from the 
simulated trajectories, as done in Ref. p^ . 

However, we can roughly estimate the number of events 
needed in the case of uniform exposure to detect the 
anisotropics predicted in our model. This can be done 
by estimating the error due to the finite number of 
events and requiring it to be smaller than cosmic variance 
and/or the deviation from predictions based on isotropy. 

Even in the case of partial sky coverage, the statisti- 
cal error on the ang ular power spectrum estimation has 
been computed |33j| . but the isotropic formula a{Ci) ~ 
Ci ~ l/(47rAfobs) {Cg is a quadratic estimator), at low £, 
is sufficient for our purposes. From Fig. 1191 we predict 
Ce ~ 0.005 above 80 EeV. Requiring cr(Cf) - 0.001 for a 
detection implies to have A'^obs ^ 80 events above 8 x 10^^ 
eV, which has not been yet reached by AGASA, but will 
be easily reached by the Pierre Auger project. 

As for the auto-correlation function Af{0), we predict 
an excess A/" ~ 4 for < 0o = 2° at 80 EeV, see Fig. [TTI 
For A^obs events, the Poissonian error on Af is cr(A/') ~ 
(^obs^o/4)-i/2. Requiring a{JV) < 0.8 for a detection 
leads to A'obs ^ 70 events needed above 8 x 10^^ eV. 

Overall, the anisotropics arising in the framework of 
this model are largely detectable by the Pierre Auger 
Observatory. 

VI. CONCLUSIONS 

We performed numerical simulations of ultrahigh en- 
ergy nuclei propagation in a structured, magnetized uni- 
verse. The sources, of density ~ lO'^Mpc"^, were dis- 
tributed proportionally to the baryon density obtained 
from the large scale structure simulation of Ref . We 
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considered either negligible magnetic fields or fields ob- 
tained from that same simulation of large scale structure 
formation and seeded by the Biermann battery around 
cosmic shocks. Their strength has been normalized to 
~ /xG in the core of the most prominent simulated col- 
lapsed structure We injected iron and protons at the 
sources, and in the case of iron we followed all secondary 
nuclei produced by photo-disintegration. 

Below the GZK feature, protons tend to dominate the 
spectrum even for exclusive iron injection, due to photo- 
disintegration of high-energy heavy nuclei. However, the 
predicted spectrum does not match the observations by 
AGASA and HiRes at these energies, and a light compo- 
nent from injection of light nuclei or protons must con- 
tribute to the spectrum at sub-GZK energies. 

Above ~ 30 EeV existing data are poor and consistent 
with an UHECR flux exclusively resulting from iron in- 
jection. As the maximum energy reachable with shock 
acceleration increases with the charge of the nucleus, it 
makes sense to expect a heavy nuclei component above 
the GZK energy. Within our scenario we found the fol- 
lowing generic predictions in the case of pure iron injec- 
tion: 

• The average mass (A) of observed particles strongly 
depends on the scenario and typically is 20 — 30 
above ~ 100 EeV. If (A) > 35 is observed at ener- 
gies below 30 EeV, the effect of extra-galactic mag- 
netic fields on propagation cannot be too strong 
since it would increase photo-disintegration. Below 
~ 100 EeV, (A) increases with increasingly steeper 
injection spectrum. 

• If some sources are located within a few Mpc from 
the observer, a flattening of the spectrum after the 
GZK feature results. This flattening is more pro- 
nounced than in the case of proton primaries, but 
cosmic variance is much larger as well. When the 
spectrum is normalized around 30 EcV, the flatten- 
ing tends to be stronger for significant EGMF. 

• A low source density and considerable magnetic 
fields, especially around the sources, predict a sig- 
nificant clustering of UHECRs at super-GZK en- 
ergies, typically on scales of order of the angular 
size of the magnetized region around the sources. 
In our simulation it is of order ~ 10°. This can be 
seen from both auto-correlation function and angu- 
lar power spectrum measurements at various ener- 
gies. 

• The injection spectra which allow to fit the ob- 
served spectra tend to be somewhat harder than 
in the case of protons, and cover the range 2.0 ^ 



a ^ 2.4, but still consistent with expectations from 
theory of shock acceleration . 

The main feature of these simulations compared to 
previous work is the distribution of structured magnetic 
fields, compatible with existing data on extragalactic 
fields. These fields increase the deflection of nuclei, and, 
as a consequence, smooth the GZK feature in the energy 
spectrum as well as the auto-correlation signal at small 
angles. We also found that, for iron injection, there is a 
major cosmic variance at the highest energies for observ- 
ables such as energy spectrum, composition and angular 
distribution, mostly due to the uncertainties in the loca- 
tion of the sources. 

At the highest energies the observed flux will be dom- 
inated by only a few sources within the GZK-distance, 
~ 50 Mpc, from the observer. Although extra-galactic 
magnetic fields can considerably smear out the images 
of the sources, their location may still be determined by 
next generation experiments to within 10 — 20°. 

In the presence of EGMF auto-correlations are reduced 
by factors of order 10, and thus more strongly than in the 
case of proton injection 14]. However, even in the pres- 
ence of the relatively strong EGMF considered in the 
present work, auto-correlations remain significant. As 
a consequence, if no auto-correlation is observed, the 
source density must be larger than 10~^Mpc~^. Recall 
that this density is motivated by the claimed AGASA 
clustering IMHa]. 

We have now explored an extensive series of plausible 
scenarios in the context of an astrophysical UHECR ori- 
gin. The scenario mostly studied here, because of the 
large magnetic fields involved over extended regions, of 
the small source density compared to interaction lengths, 
and of the fact that we inject iron nuclei at the sources, 
appears to be the "worst case" in terms of fluctuations 
of the predictions of observables. Hopefully, the Pierre 
Auger Observatory will soon allow to select a speciflc 
scenario if its data is compatible with an astrophysical 
origin of UHECRs. 
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